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Abstract 

Let r be the first hitting time of the point 1 by the geometric Brownian motion 
X(t) = xexp (B(t) — 2 fit) with drift p, ^ 0 starting from x > 1. Here B(t ) is the Brownian 
motion starting from 0 with E°B 2 (t ) = 2t. We provide an integral formula for the 
density function of the stopped exponential functional A(t) = fj X 2 {t)dt and determine 
its asymptotic behaviour at infinity. Although we basically rely on methods developed 
in Ena, the present paper also covers the case of arbitrary drifts ^ ^ 0 and provides 
a significant unification and extension of results of the above-mentioned paper. As a 
corollary we provide an integral formula and give asymptotic behaviour at infinity of the 
Poisson kernel for half-spaces for Brownian motion with drift in real hyperbolic spaces of 
arbitrary dimension. 
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1 Introduction 


Let B(t) be the standard Brownian motion starting from 0 and normalized such that E°B(t) 2 = 
2 1. Consider the following linear SDE: 

dX(t ) = X(t)dB(t ) — (2/i — 1 )X(t)dt, Y(0) = x > 0, // ^ 0. 

The strong unique non-exploding solution is given by 

X(t) = xexp (B(t) — 2 fit). 

The process {X(t);t ^ 0} is called a geometric Brownian motion or exponential Brownian 
motion and along with the additive functional from the process X(t), 

A x (t) = f X 2 (t)dt= I x 2 exp 2(B(s) — 2/is)ds 
Jo Jo 

is of primary interest in mathematical finance and insurance theory (see, e.g., jDj, |GY| or IY2F . 
Also there is a connection of the above functional with the Brownian motion in hyperbolic 
halfspaces (see eg. jY2j, pUFj . jBCFYj and Section 5). 

The distribution of A x (t) for fixed 0 < t < oo has been a subject of study in substantial 
number of papers (see e.g. 0 , na, isa. [M]). In this paper we investigate the properties of 
the density function q^ of the (stopped) additive functional A x (r) where r is the first hitting 
time of the point 1 by the process X(t) (starting from x > 1). From the strong Markov 
property it easily follows that the distribution of A x (t) is closely related to the distribution 
of Ae(oo). The latter one is astonishingly simple; it is identical with the distribution of the 
random variable x 2 /4Z , where Z is a T(/i, l)-distributed random variable (with the density 
r(/i) _1 tG _1 e _ “). This fundamental result is due to D. Dufresne [13] (see also |Ylj . ;Y2]) and is 
of primary importance here. Knowing the Laplace transform of Aj,(oo) one can easily derive 
(via the strong Markov property) the form of the Laplace transform of the random variable A(t) 
as a suitable ratio of Bessel functions. For the reader’s convenience we present this argument 
in Preliminaries. 

On the other hand the distributions of A x (r) and A x (oo) are closely related to hitting times 
of Bessel processes. The main fact here is Lamperti’s representation which reads that there 
exists a Bessel process w ith index —/x, starting from x > 1, such that the process X(t) 

admits the following representation (see JI] and Ex. 11.1.28 in [KYI ): 

X(t) = (A x (t)), t> 0. 

We refer the reader to IRYl for the account on Bessel processes (see also (sni). From the 
Lamperti representation it follows immediately that A x (r) and A x (oo) can be regarded as 
hitting times of 1 and 0, respectively, of the Bessel process starting from x > 1. Also it is 

possible to relate A x {r) and Ae(oo) to last exit times for appropriate Bessel processes. There is 
a vast existing literature on that subject. Such hitting or last exit times were studied in papers 
by Getoor |H] . Getoor-Sharpe jlGSI) . Kent jKj and Pitman-Yor PY . In Getoor [G], Kent [K] 
and Getoor-Sharp [GS| the Laplace transform of the distribution of A{r) is derived as a ratio 
of Bessel functions (see Preliminaries). For an exhausting discussion on that subject we refer 
to HUYj. 

Our main focus in this paper is to provide an integral formula for the density of A x (t) (see 
Theorem ESI)- In the case p — 1/2 this density is well known to be 1/2-stable subordinator. 
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Kent in jK] writes that it is possible to obtain the explicit form of the density for p = 3/2 but 
he does not provide any formula nor details on that. 

At this point let us mention that various ratios of Bessel functions were proved to be 
completely monotonic functions hence they are Laplace transforms of probability distributions 
(see Ismail jHj,jEl, Ismail-Kelker OKI . Kent lK|). For a survey on this theme see Pitman and 
Yor and also dt Regarding the case which is considered in our paper Ismail and Kelker 
ITKj showed by purely analytical methods that it is an infinitely divisible distribution. 

The main purpose of the paper is to obtain a suitable representation for the density function 
of the functional A x (t) along with its asymptotic properties, for arbitrary drift terms p ^ 0. 
Then we apply these results to derive an integral representation of the Poisson kernel for 
subspaces on real hyperbolic spaces, for hyperbolic Brownian motion with arbitrary drift, which 
extends and simplifies the results and proofs from [BGSj . 

The paper is organized as follows. In Preliminaries we collect basic information needed in 
the sequel. 

In Section 3 we obtain a suitable representation of the density of the functional A x (t). For 
this purpose we apply an integral representation given in | IBGS| for ratio of Bessel functions. 

In Section 4 we exhibit the exact asymptotics of the density of A x (t) at infinity for all 
drifts fi ^ 0. Again, we essentially follow the idea of CUE!. However, applying more direct 
probabilistic arguments, we are able to simplify our presentation considerably. 

In Section 5 we show how to apply results obtained in preceding sections to obtain repre¬ 
sentation and asymptotic properties of Poisson kernel of subspaces on real hyperbolic spaces 
of arbitrary dimension, for Brownian motion with drift. Some of these asymptotic properties 
were studied in in the case of dimension 2, where it was shown that the distribution 

of Poisson kernel belongs to the stable domain of attraction. Our asymptotic results may be 
viewed as the extension of those obtained in ihiif!. For a related result see also jBGFYj . 


2 Preliminaries 

Let 0 < a < x and let r a be the first, hitting time of the point a by the geometric Brownian 
motion with drift fi Y 0 starting at x: 

r a = inf{Y > 0; xexp(F>(t) — 2/it) = a}. 

The fact that r < oo a.e. follows from the property inf i>0 R(t) = — oo. 

Further, define 

A x (t) = / exp 2(£>(t) — 2fit)dt. 

Jo 

By the strong Markov property of Brownian motion we obtain 
Basic relationship (for /i > 0). 

Observe that x 2 exp 2(IF(r a ) — 2/ir a ) = a 2 , hence 

r-T a pOO 

A^oo) = x 2 exp 2{W(s) — 2/is)ds + x 2 / exp 2{W(s) — 2/is)ds 

Jo j T a 

pOO 

= A x (Ta) + x 2 exp 2(W(r a ) — 2/iT a ) / exp 2{W{s + r a ) - W{r a ) - 2/is)ds 

Jo 

= A x (T a )+A' a ( oo), (1) 

where A^oo) is a copy of A a (oo), independent from A x {r a ). 
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Dufresne [D] (see also Getoor 0 , Kent jKj, Getoor-Sharpe 1GB] . where the result is given 
in the context of Bessel processes) showed that the density and the Laplace transform of Ai(oo) 
have the following form: 


2~ 2 m e -i/4t 

= i^) 


( 2 ) 


£{^( 00 )}^) = £{l/4Z}(r 2 ) = 2 Y^K,(r), p > 0. 

From this and elementary properties of Laplace transform one immediately obtains 


( 3 ) 


C{Ax{Ta)}^) 


(xy /y(xr) 

\aJ K M (ar)’ 


( 4 ) 


where K u is the modihed Bessel function of second type with index /u. The above Laplace 
transform appears in many papers (see eg. Kent [Kj, Getoor-Sharpe ESI, Ismail-Kelker |IK| ) 
and here we followed a derivation from lEE As mentioned in the Introduction we present it 
for the readers convenience. 

From the continuity of both sides of 0 with respect to /i ^ 0 we obtain that the above 
formula remains valid also for ft — 0. 

Let us remark that for /1 = 0 the right-hand side of 0 gives also the classical formula (|Hj) 
for the Laplace transform of Brownian motion hitting time T a of the centered circle with radius 
a in IR 2 from the point i/6l 2 such that \y\ — x > a. 

By the formula 0 follows directly that for t > 0 we have 


Ata,tx{j~ta) t A ax (r a '), 


where the symbol = denotes the equality of distribution. 

Therefore, from now on we may and do assume that a = 1 and x > 1 is fixed. We write r 
instead of r\ and A(r) instead of Ai tX {ri). 

We conclude this section with the following technical lemma 


Lemma 2.1. 

poo \2s+2 

/ e~ y2/4t h s (t/X 2 )dt = 

Jo (a 2 + r 

Proof. Indeed, we obtain 


r e -.v 4. e -AV4,yy^ = 

Jo t 1+s 

2 2 s^2s+2 poo 

= (FT r{s)KMdu 


v*+ 2 r dt/tf + a 2 ) 

(V + V 2 Y Jo (f/fo 2 + A 2 ))'+* 

_ 2 2 T(s)A 2s+2 
(A 2 + y 2 ) s 


( 5 ) 


□ 

Throughout the paper we use the following convention: by c, C we always denote nonnega¬ 
tive constants which may depend on other constant parameters only. The value of c or C may 
change from line to line in a chain of estimates. 

The notion p{u) ~ q(u), u —> Uq means that the ratio q/p —> 1 when u —> Uq. 
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3 Representation of density of A(r) 

In this section we give a representation formula for the density function of functional A(t), for 
arbitrary /i ^ 0. 

From now on we use the following notation, partially introduced in the preceding Section: 


x > 1, a — 1, A = x — 1. 


and we denote by qy the density function of the functional A{r). 

We begin with stating a more general version of a lemma, taken from 
identical and is omitted. 


BGSI . The proof is 


Lemma 3.1. Let ji ^ 0. Suppose that 


X 


Q( z ) = z - (f ~ V 4 )^’ zeC. 
Define F\(z) by the following formula: 

ze Xz x^K^{xz) — x ll ~^Q(z)K fl (z) 


X F x (z) = 


KJz) 


Then 


F\(z) = 0(z 1 ), z -»■ oo. 
and there exists a function w\ such that 

poo 

F\(z ) = / e~ zv w x (v)dv. 


Moreover, 


and, for /i > 1/2, 


aF 1//2 (/i 2 — l/4)/2x = / w\(y)dv , 


2aF 1/<2 = / KW\(v)dv, 


where k = (A + u) 2 — A 2 = u(2A + u). 
For /i = 1/2 we have F\(z) = 0. 


( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 


An explicit formula for the function w\ will be provided in the sequel. 
The following formula is crucial for our considerations: 


Lemma 3.2. Let p : [0, oo) 
the following property 


M + belong to F 1 (0,oo) such that its Laplace transform (p has 


t l t 2 <p(t) G L 1 (0, oo). 


Then we have 



( 11 ) 
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Proof. Let A be a random variable with absolutely continuous distribution with the density 
function <p and let W(t) be Brownian Motion starting from 0 (such that EW 2 (t) = 2 1), inde¬ 
pendent from A. It is easy to see that the value of the Laplace transform of A at the point r 2 , 
that is (p{r 2 ) is equal to the Fourier (or cosine transform) of W(A) at the point r. Observe that 
our assumption assures that this Fourier transform belongs to L 1 (0, oo). Thus, the left-hand 
side of m is the inversion formula applied for the Fourier transform and gives the density 
of W(A). The right-hand side results from the direct computation of this density, taking into 
account the independence of W and A and the particular form of the (gaussian) density of 
W(t). ' □ 

We are ready to state our representation formula. 

Theorem 3.3. For p ^ 0 we have 

p —X 2 /At f ro o \ 

q^t) = A-^=r- ( x^ 1/2 /2t + J (e~ K ' At - 1) w x (v)dv) . (12) 

For p > 1/2 we have 

-\ 2 /4t roo 

quit) = A— -j=- / (e _K/ * - 1 + k/A t)w x (y)dv. (13) 

VTTt Jo 

Proof. The proof relies on application of the formula m for the function 


( p(r 2 ) = x ^ 


K^,(xr) 

K»(r) 


Note that by © and 0 it follows that 
.K^xr) 


rx r 


K,{r) 


= e~ Xr XF\(r) + e~ Xr x»- 1/2 Q(r ) 

= e~ Xr x»- ll2 r + e" Ar AF A (r) - e~ Xr x^ 12 ^ 2 - l/4)A/2x 


= e~ Xr x»- ll2 r + Ae _Ar / (e 


l) w\(v)dv. 


We recall that Q(r) = r — (/i 2 — l/4)A/2x and F x (r) = J 0 °° e rv w x (y)dv. To simplify the proof 
we introduce new notation. Namely, denote 


wf(v) = - 


dwf(v ) 
dv 


w x (v)dv , that is 
Then, by integration by parts we obtain 


= w x (v). 


(e rv — l) w x (v)dv = (e rv — l) w* (n)|^° + r / e rv wf{y)d 


= r 


e rv wf(v)dv. 
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Thus, the left-hand side of the formula CP is now of the form 

t K n(xr) 


rx r 


■ cos (ry)/rdr 


Jo K„(r) 

poo 

x^~ 1/2 / e~ Xr cos(ry)dr 


'0 


r* oo 


"OO 


+ / Ae Ar { / (e — l) w\(v)dv} cos(ry)/rdr 

Jo Jo 

poo 

= x / e” Ar cos (ry)dr 


r»oo /»oo 

^-Xr 


+ / Ae Ar { / e TV wf(v)dv}cos(ry)dr 

Jo Jo 

Using the standard formula for the Laplace transform of the cosine function we obtain 

A 


4 1J (!/)=^ 1/2 


and 


A 2 +1/ 2 ’ 


.do (A + v) + 1/ 


By the form of dd^(y), © (applied for s = 1), © and (HU we obtain 


H£\y) = 


( 2 ) 

Analogously, for H^ ' we obtain 

A / ‘ 0 ° 


sffe) = 


2 2 Jo 

A /■“ 


e ~y 2 / 4t I / (A + u)wf(u)e- (A+ " )2/4t du 


'0 

/•oo 


dt 

¥ 




Observe now that 

d(- e -( A +") 2 / 4 *) # 


dw 


wf(u)du = -e“ (A+,;)2/4t u;f(u)|“+ / e - (A+,;)2/4t 


dwf(i 


dv 


-dv 


= e~ x2/4t wf( 0)+ / e~ < ' X+v ' )2 ^ 4t w\(v)dv. 


We also have 


»f( 0) =- r^= 


du 


w\(v)dv. 


(2) 

By the above identities and the form of dd^' we obtain 

A ro ° 


^i 2 ) (d) = 77 / e y2/4t / ( e K/4t - l)w A (n)du —. 


dt 
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Combining the above identities we obtain that the left-hand side of m takes the form 

H^y) = e-v 2 / At e- x2 ' At (x^ 2 ,/2t + J"~ lH (v)di^j j. 

Taking into account the right-hand side of m and continuity of q M with respect to t (see 
properties of w\ below), we obtain (fT2j) . When /i > 1/2, then using (fTUl) we obtain m- 

□ 


Below we give a description of the function w\. We rely here on results contained in jBGSj . 
The formulas depend on the zeros of the function K^(z). 

Even if in general the values of these zeros are not given explicitly, we are able to prove 
some important properties (as boundedeness or asymptotics) of w x , which are essential in 
applications. Moreover, for some values of /i we provide explicit formulas as well (see Corollary 
13.51) . 

The function K tl (z) extends to an entire function when p — 1/2 is an integer and has a 
holomorphic extension to C\ (—oo, 0] when p — 1/2 is not an integer. Denote the set of zeros of 
the function K tJ (z) by Z — {zi, ..., Zk tl }• We give some information about these zeros (cf. [Ej, 
p. 62) needed in the sequel. Recall that k^ = ji — 1/2 when p — 1/2 e N. For p — 1/2 ^ N, k M 
is the even number closest to p — 1/2. In particular, for 0 ^ p < 3/2 we have k^ = 0; for ft = 2 
and 3 we have k^ = 2. The functions K M and K /x _ 1 have no common zeros. 

As in urns], we need an additional notation to describe the function w\. Define for /i > 0 




yk Zie Xzi K^xZi) 2] 
A “ K^^Zi) 


(14) 


When p + 1/2 </ N and p > 0 we define 


W 2 ,\(v) = -cos(7r/i) 


x H 


Ifi {xu) K^u) - I^{u)K^ (xu) 


A J o cos 2 (7r/i)/W(w) + (7 tI^(u) + sin (7 x^K^u)) 


-e~ Xu e~ vu udu. 


(15) 


We now formulate our representation theorem for the function w\. The proof of the main 
part is the same as in ]IBGS| and is omitted; we only show asymptotic properties of the function 
w 2 ,\- For fi = 0 this is new; behaviour for fi = (n — l)/2 was shown in jBGSj . Nevertheless, we 
present here a new and unified proof based on tauberian theorems. 


Theorem 3.4. In the case fi — 1/2 e N 


w\(v) = w 1>x (v); 


while, in the case when fi — 1/2 ^ N 

w\{v) = w ljX (v) + w 2 ,x(v). 

Moreover, we have sup,,^ 0 |iu*(i;)| < 00 , and 

— cos(tt/j,)w 2 ,\(v) ^ 0, v ^ 0, (/i — 1/2 N); 


v k \wi : \(v)\dv < 00 , k — 1,2,...; 

) 

lim v k Wi )X {y) = 0, A: = 1,2,...; 











, o.io . . — COs(7Tu)r(2u + 2) x 2fl — 1 . , , . 

1m, „ w,Av) = 2 4i44 + 1) A ■ - 1/2 * N ' " > 0); 


lim ( vlogv) 2 W 2 t \(v ) =-- log a;, /or fi = 0 

»oo ’ 2A 


Proof. Denote 


4,a( m ) = 


_4 (to) A^(m) - I^(u)K^ (xu) 

cos 2 44A" 2 (4 + 44(4 + sin 44 ^4 (4) 


-e~ Xu u. 


( 16 ) 


Observe that the function 4 ,a (4 is non-negative (I^(u)/K^(u) is increasing for u > 0) and 
w 2 ,a(v) is the Laplace transform of — cos44x4,a at v. We claim that for u —> 0+ we have 
the following asymptotics for 4 ,a' 


KA U ) 


0+4(1 — x 2fJ ’)u 2fJ,+1 , 

C M 

u logx(logn) -2 , 


for /i > 0, 
for p = 0. 


(17) 


Applying Karamata’s Tauberian theorem (see, e.g. jFe|) we obtain that the asymptotic be¬ 
haviour of W2,\( v ) i s 


lim w 2 ,x ( v ) 

KX) 



— C0S(7T/i) 

A 


r(2 + 2/i)x M . 


This together with (El implies: 

lim v 2 ^ l+2 W 2 ,\(v) = 

v —>00 

lim (v logv) 2 w 2 A v ) 
^—>■00 


— COs(7T/i)r(2p + 2) o4 — 1 

2 2 ^- 3 r(/i)r(/i + 1 ) a 

= _ 2A logX ’ /i = 0 ' 


4 -l/2^N, fi> 0), 


To prove dH we apply the following asymptotics. When u —> 0 we have: 

4(“) ~ c +4 ^4(4 ~ 

4(4 = l + o(l), K 0 (u) = log(2/44(4 + -0(1) + o(l), 

with = 2 -/ yr(p + 1)) 4 = 2 M_1 T(4 and where if is the Euler function. Then 


(18) 

(19) 




4, a 4) = 


4 (xu) i4(4 


AT 2 4) + (7r4(w)) 2 + 27rsin(7r/i)J4(M)4(u) 


4 ( CT ) 


x _ 4 (u) 14 (to) 

4 4414 (4 


4+z+4(i — x 2/i )> f° r + > 0 , 

M for 11 — 0. 


(log u) ; 


2 j 


444444 \ 

4 44 K^(u)) 


□ 


Examples. 

To illustrate representation theory developed so far we write down explicit integral formulas 
for the density q /t in some special cases of /j. All formulas appearing here follow directly from 
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Theorems 13.31 and ro If0^/r<3/2 then w x = in 2 ,A and functions w x have constant sign. For 
/i = 0 and /i = 1 the function w x has simpler form, which we exhibit here. If /i = 1/2 then the 
function Fa = 0 and q,,{t) reduces to the standard 1/2-stable subordinator. For /i + 1/2 e N 
we have, in turn that w x = W\. x and the form of W\ t \ can be computed calculating residues of 
simple rational functions (see the formula for w\ jX or calculations in |BGSj ). Again, we write 
the explicit form of w x for /i = 3/2 and /i = 5/2. 


Corollary 3.5. If fi = 0 then 
-w x (v) = 

and 


1 f°° I 0 (xu)K 0 (u ) - K 0 (xu)I 0 (u) ^ 


Ad o Ko{u) + n 2 Io{u) 


-e- uA e~ vu udu, 


3 -a 2 /4 t 


y/irt 


%{t) = A 7 = (A + 1) 1/2 / 2i + / (1 - e K/4i ) (-w A (r))(i? 


If H — 1/2 i/ien 


?/*(*) = A 


e~ A2 / 4t 
2 a/ 7rf 3 


If /a = 1 t/ien 


wa(w) = 


A + 1 Z 00 /iOro)Ad(M) - K\{xu)Ii(vl) _ uX _ 


A 


Kl(u) + 7T 2 /{(m) 


e- UA e~ vu udu, 


and 


-,—\ 2 /4t roo 


qn(t) = A ' — / (e K/4i - 1 + K/4f)u;A(w)dn. 
Vvrf Jo 


If li = 3/2 t/ien w x (v) = e v and 

D —\ 2 /4t roo 


q^(t) = A —-=^ / (e K,/4f — 1 + /c/4 t)e V dv. 
Vvr t Jo 


If fi = 5/2 i/ien 


w x (v) = 3e 3 w// 2 [(2A + 1) cos(a/3w/ 2) + \/3sin(\/3i;/2)], 


and 


q^(t) = A 




1 + k/A t)w x {v)dv. 


4 Asymptotic behaviour of A(r) 

In this section we prove the following 
Theorem 4.1. The density of A(r) satisfies: 

lim F + V(t) = Cp, if li >0, 

t—>oo 

lim (logt) 2 %(t) = C 0 , if fi = 0 , 

t —»00 

for some positive C M . 
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Lemma 4.2. Let /a > 0. There exists a constant C > 1 such that 

C~ l r» < P(A(r) >t)< Ct~t > 1. 

Proof. From © we infer that (x > 1, a = 1): 

P(A( oo) > t) < P(A(t) > (1 — x~ l )tx 2 ) + P(A( oo) > xt ), 

which implies 


P(t < H(oo) < xt) < P(A(t) > (x — l)xt). (20) 

Moreover, 

P(A(t) > x 2 t ) < P(H(oo) > £). (21) 

Now the lemma follows from (I2U1) and (I2T1) since P(A(oo) > t) ~ ct~ M for some c > 0 by (J2J). □ 

Recall that k = (A + v ) 2 — A 2 . 

Lemma 4.3. Let m e N fee such that 2<m</i+l/2. T/ien 

/ °° i / k \ j i—Il m t 100 

P) )*- = / K m Mv)dv = C m . (22) 

Moreover, C m = 0, for 2 < m < /i + 1/2, and, in the case m = fi + 1/2 e N, we /lave C m > 0. 
Proof. Denote 

By elementary calculations 


3 -K/4t 


1 


E (-Dh 




K, \ J 

4t/ 


< 


m! 


At, 


Hence 

K m 

\^{t,m)\ < |w A (n)| —. 

ml 

Under the assumption on m the function \w\(v)\n m is integrable so the formula (1221 ) follows 
from the bounded convergence theorem. 

Suppose that there exists C m 0, with m having properties as above, and denote 
m 0 = inf{m eN: 2<m^/i+l/2, C m ^ 0}. 

Then from the first part of the proof we have 

lim t mo+1/2 q^{t) = C mo > 0, 

t—>oo 


which implies that 

lim t m °~ 1/2 P(A(r) > t) = C mo (m 0 - 1/2) -1 . 

t — >oo 

From Lemma n~ 2 i we infer that C mo > 0 if and only if mo — 1/2 = fi. In particular, we then 
have /i + 1/2 e N. We also obtained that ifm</i + l/2 then C m = 0. This completes the 
proof of the lemma. □ 
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Remark. The above lemma yields, in particular, that for m e N, 2^m</x+l/2 

/»oo 

/ K m w\(v)dv = 0. 

Jo 

Thus, the representation formula for the density can be written for /x ^ 1/2 as follows: 

„-A 2 /4t roc 


Qii(t) = A- 




(e 


K/4t - Z (- 1 )' 




( 23 ) 


where l = [/x + 1/2], if p — 1/2 ^ N, and l — /x — 1/2 otherwise. 

We now prove our theorem. 

Proof of Theorem mu 

Proof. For /x = 1/2 the density has a particularly simple form (see Corollary 13.51) and the 
theorem clearly holds true. Hence, we assume throughout the remainder of the proof that 
/x ^ 1/2. Next, if /x — 1/2 6 N then the Remark above together with Lemma 14.31 yield our 
theorem at once. 

Thus, we assume for what follows that p — 1/2 ^ N and let l = [/x + 1/2]. Denote 

" - V (-l)T (£) V- 


m = 


o sSjsSZ 


We prove that 
and 


lim I(t) — C > 0, if /x > 0; 


£—>•00 


lim (logt) 2 /(t) = C > 0, if fi = 0. 


£—>oo 


(24) 

(25) 


Applying change of variable k = 4sf we obtain 

v = \/4st + A 2 — A, dv = 


2 1 


y/Ast + A 2 


ds, 


so 


£(*) = / ^x(s,t)(e~ 


E (-i)A)^. 

oyjyi • y ' 


where 


V’a(S) 0 = + P - A)- 


2t /x+3 /2 


4 / 4sf + A 2 

We claim that for /x > 0 and t, s G f? + there is a constant C > 0 such that 

IVaM)I ^Cs-^ 2 . 

For /x = 0 and t, s G f? + our claim is: 

(logtfliMs,f) < Cmax{l, (logs) 2 }s _3/2 , 


(26) 


(27) 


where C > 0. The above claims prove the relations (l^H) and (1251) . Indeed, consider first the 
case /x > 0. Then the absolute value of the expression under the integral I(t) can be estimated 
by the integrable function 


c 


E 

o^Ki J ' 


s -v-3/2 < c nrinjl, s}s z s _A1_3/2 = Cminjl, 
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where 5 = /a + 1/2 — [/a+1/2], 0 < S < 1. Then the proof of (1%H) is concluded by passing 
t —> oo and using the asymptotics of w\ (see Theorem 13.41) . 

Now, consider the case n — 0. Observe that the absolute value of the integrand in (log t) 2 I{t) 
is estimated by the integrable function 


Omax{l, (logs) 2 }s 3 ^ 2 |e s — 1|, so 

lirn (log t] 2 I(t) = f (l - e~ s ) s~ 3/2 ds. 

t_>0 ° 2.X J 0 

Here we take into account 


(log t) 2 -if x (s,t) « 


log a; (logf) 2 log a; 

2A s 3 / 2 (logsf) 2 2As 3 / 2 ’ 


oo, 


by using the asymptotics of w\ (see Theorem 13.41) . 

Now we can conclude the proof of the asymptotic behaviour of q^ft) in the case /a —1/2 ^ N. 
Note that for 0 < /a < 1/2 we have 


= A 


= A 


e -A 2 /4 1 

\[tP 

g -A 2 /4 1 


7r 


(tV~ 1/2 /2 + ++> jf (e~ K/4t - 1) wxWdv'j 
(^- i / 2 x m - i/2 / 2 + _+ xc /^ - 


where C is the constant from the right-hand side of the formula (1^1) . The same argument for 
/a = 0 shows that 

(log tftq^(t) -> \Cfy/n, 

where C is the constant from the right-hand side of the formula m- For /i > 1/2 the 
asymptotics of q tl (t) directly follows from the formula (|24|) . □ 

We now justify our claims (OBI) and (EH)- We use notation as introduced in the proof of the 
theorem. 


Lemma 4.4. If /i> 0 there is c> 0 such that 

MM)| ^ cs~»- 3/2 . 

For ja —— 0 there is C > 0 such that 

(logt) 2 |'0 A (s,t)| < C max{l, (log s) 2 }s _3/2 . 

Proof. We begin with EHD first. Let /i > 0. Using boundedness of the function w\ we can 
estimate for ts < 1: 

3/2 

MM)| < sup \w\(v)\—=== < Cs~^ 3/2 . 
v>0 V4 st + A 2 

For ts > 1 we use the asymptotics of w\ at oo (see Theorem EH) to arrive at 

MM) I < C^ +3 / 2 ( S f)-^ 3/2 = Cs-^~ 3/2 . 

For /a = 0 we estimate our expression as follows. Again using boundedness of the function ui\ 
we estimate for ts < 2: 

a 3 / 2 

\fj x (s,t)\ < sup\wx(v)\^=== <Cs 3/2 . 
u>o V4 st + A 2 
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For 2 < ts we use the asymptotics of w\ to get 

»••»•'**»'«»R?* 0 ' iv'j. A 

Next observe that for 2 < ts < y/t we have (logf) 2 < 4(logs) 2 . 


□ 


5 Hiperbolic Brownian motion with drift 

Consider the half-space model of the n-dimensional real hyperbolic space 

HP = {(xi,..., x n _i, x n ) G IP” 1 xl: x n > 0}. 

The Riemannian metric, the volume element and the Laplace-Beltrami operator are given by 

ds 2 = 


2 dx\ + ... + dx 2 _ x + dx 2 


X r 


dV = 


dxi...dx n -idx r 


rpT i 

•hr 


A = x2 nC£2 - ( n - 2)x n 0 n , 


i =1 


respectively (here di = i — 1,... ,n ). For // > 0 let a = 2/i — n + 1. We also introduce the 
operator: 

n 

= A - ax n 5 n = x 2 (^ Of) - (2 /i - l)x n d n . 


i=l 


Let {Bi(t))i=i... n be a family of independent classical Brownian motions on M with the generator 
-d ~2 (and not i.e. the variance E°Bf(t) = 2 1. Then the Brownian motion on HP, X = 

(Xi) i= i_ n can be described by the following system of stochastic differential equations 


dX^t) = 
dX 2 (t) = 


X n {t)dBi{t) 

X n (t)dB 2 (t) 


dX n (t) = X n (t)dB n (t ) - in - 2 )X n (t)dt. 

More generally, if we replace n — 2 by a + n - 2 = 2/i — 1 that the corresponding process will 
be called the Brownian motion on HP with drift a. By the Ito formula one verifies that the 
generator of the solution of this system is A /t . Moreover, it can be easily checked that the 
solution is given by 

X^t) = W(0) + X n (t)dB 1 (s) 

X 2 (t) = X 2 (0 ) + J*X n (t)dB 2 (s) 

x n (t) = X n (0) exp(B n (t) - 2/xt). 

Dehne the projection': R n 3 u = (ui,...,u n ) —> u = (u±, ..., u n _i) G HP -1 . In particular, 
X(t) = (Ai(t),..., A n _i(t)). From the representation above one may easily verify ( e.g. by 
computing Fourier transforms) that 

X(t) = A(0) + B(f Xl(s)ds), t > 0. (28) 

Jo 
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where B(t) = (i?i(f),..., B n _i(t)) is n — 1 dimensional Brownian motion independent of the 
process X n (t). Consider a half-space D = {u G HP : u n > a} for some fixed a > 0. To simplify 
the notation we choose a = 1. Define 


r = inf{i ^ 0 : X(t) ^ -D} = inf{t ^ 0 : X n (t ) = 1}. 


By Pi(u,y), u = (mi, w 2 , ■■■,u n ) e D, y = (■ ?/i,t/ 2 , ■■■iVn-ii 1) £ dD we denote the Poisson kernel 
of D, ie. the density of the distribution of X(r) starting at u (since X n (r) = 1 it is enough to 
consider the distribution of X(r)). From (f28l) it is obvious that 


X(t)^u + B(A(t)), 

where the functional A(r) (starting from u n > 1) is independent of B(t). For further con¬ 
siderations we may take u = 0 and u n = x > 1 , so the starting point of the n dimensional 
process X(-) is (0,..., 0, x) G D. Since A{r) and B(t ) are independent, we have the following 
representation of the Poisson kernel: 


Corollary 5.1. 

V) = (i 7 !) l-rn l 

Observe that for p — 1/2 the functional A(r) has the standard asymmetric 1/2-stable 
distribution and the Poisson kernel is exactly n — 1 dimensional Cauchy density, so in what 
follows, we consider p ^ 1/2. 

For the sake of simplicity we write \y\ as p. Taking into account the formula (( 21 ) and m 
we obtain the following representation of the Poisson kernel of the set D 

Theorem 5.2. For 0 ^ p < 1/2 we have 




(n — 2)(A + l)^ -1 / 2 — 


w\(v)L#(\, p, v)dv 


27 T n / 2 (A 2 + p 2 ) n / 2 
where A, p, v ) is the following function 

L#( A, p, v ) = (A 2 + p 2 )[((A + v) 2 + p 2 )^ 2 " 1 - (A 2 + p 2 ) n/2 “ 1 ]. 
For p > 1/2 we obtain 

p , X r (f ~!) A r°° w x (v)L(X 1 p,v)dv 

1[X,V) 2vr-/ 2 (A 2 + p 2 ) n / 2 J 0 ((A + n) 2 + p 2 )t- 1 ’ 


1 0 ((A + n) 2 + p 2 )t- 1 


with L(X,p,v ) defined by 


L(X,p,v ) 

= (n/2 - 1)((A + u) 2 - A 2 )((A + n) 2 + p 2 )^ 2 " 1 - (A 2 + p 2 )[((A + n) 2 + p 2 )^ 2 ” 1 - (A 2 + p 2 )^ 2 ” 1 ]. 
Proof. By the formula © we obtain 

r = 22,r < s ) 

Jo ( 1+ * (A 2 + p 2 )*' 

Applying the above formula (with different constants) for three terms appearing in the repre¬ 
sentation of q M we obtain, 
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and 


and 



e -p 2 /4tg-A 2 /4t 


dt 

■f-l+n/2 


2 r T(f) 

( A 2 + p 2)n/2’ 


e -P 2 /4t e -A 2/ 4t 


'0 


r( i - J > 2 


n —2 


e " /4t iu A (i;) 


w\{v)dv 


dt 
t n ! 2 


/o ((A + n) 2 + p 2 )t 


,2\§-l ’ 


r»00 COO 

g-p 2 /4t g -A 2 /4t / 


W A (w) 


dt 

W 2 


r /™ _ non-2 _Wx(jj)d^_ 

2 j y 0 (A2 + P 2)f-i 


Taking into account appropriate constants, we obtain the formulas for the Poisson kernel. □ 


Theorem 5.3. 

lim | 2 /| n+2/i_1 Pi(a:,y) = c M , p > 0; 

|y|—>oo 

lim (log \ y\) 2 \y\ n ~ 1 Pi(x, y) = c 0 , p = 0; 

!sh<» 

for some positive c M . 

Proof. From Corollary 15.11 Theorem 14.11 and arguments presented in its proof (boundedeness 
of t^ +1 q^{t) for /i > 0 or (log t) 2 q^(t), for /i = 0) we obtain that the asymptotic behaviour of 
Pi(x, y) is the same (up to a positive constant) as of 



£-l-p-(n-l)/2 e -|y| 2 /4t^. _ 22 p+n-l | ^ |-2p-n+l 


M P+(n-3)/2 e -u du ^ 


for p > 0. Similar arguments apply for p — 0. 


□ 
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